Dynamics of Black Hole Pairs I: 
Periodic Tables 
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Although the orbits of comparable mass, spinning black holes seem to defy simple decoding, 
we find a means to decipher all such orbits. The dynamics is complicated by extreme perihelion 
precession compounded by spin-induced precession. We are able to quantitatively define and describe 
the fully three dimensional motion of comparable mass binaries with one black hole spinning and 
expose an underlying simplicity. To do so, we untangle the dynamics by capturing the motion 
in the orbital plane and explicitly separate out the precession of the plane itself. Our system 
is defined by the 3PN Hamiltonian plus spin-orbit coupling for one spinning black hole with a 
non-spinning companion. Our results are twofold: (1) We derive highly simplified equations of 
motion in a non-orthogonal orbital basis, and (2) we define a complete taxonomy for fully three- 
dimensional orbits. More than just a naming system, the taxonomy provides unambiguous and 
quantitative descriptions of the orbits, including a determination of the zoom-whirliness of any 
given orbit. Through a correspondence with the rationals, we are able to show that zoom- whirl 
behavior is prevalent in comparable mass binaries in the strong- field regime, as it is for extreme- 
mass-ratio binaries in the strong- field. A first significant conclusion that can be drawn from this 
analysis is that all generic orbits in the final stages of inspiral under gravitational radiation losses 
are characterized by precessing clovers with few leaves and that no orbit will behave like the tightly 
precessing ellipse of Mercury. The gravitational waveform produced by these low-leaf clovers will 
reflect the natural harmonics of the orbital basis - harmonics that, importantly, depend only on 
radius. The significance for gravitational wave astronomy will depend on the number of windings 
the pair executes in the strong-field regime and could be more conspicuous for intermediate mass 
pairs than for stellar mass pairs. The 3PN system studied provides an example of a general method 
that can be applied to any effective description of black hole pairs. 



At first glance, the orbits of a black hole pair resist 
coherent description. Dynamically, black hole pairs in- 
volve non-linear relativistic effects leading to an extreme 
form of perihelion precession, coined zoom-whirl behav- 
ior, as well as spin precession that in turn drives orbits 
out of a plane. The three-dimensional precessions fill 
out a tangled path that shapes the gravitational waves 
both LIGOjiLand LISA || were designed to observe 
@, i, Si, 011, S- Despite appearances, we show the 
path can be untangled and a coherent description of fully 
three-dimensional precessing orbits proves entirely possi- 
ble for one spinning black hole, one non-spinning. 

The spacetime around two orbiting black holes eludes 
analytic description. While the impressive breakthroughs 
of numerical relativity [TO, [O, E, El, Q E H promise 
to resolve the final plunge of a black hole pair, compu- 
tational expense relegates the majority of the inspiral 
to analytic methods. Several groups have gone to great 
pains to build a Hamiltonian formulation for two black 
holes in a Post-Newtonian (PN) approximation for ana- 
lytic computation of orbits and the gravitational waves 
they generate. As a contribution to this great campaign, 
we use the conservative 3PN Hamiltonian plus spin-orbit 
couplings to deconstruct the full three-dimensional dy- 
namics of black hole pairs, in the absence of radiation re- 
action. Our binaries are composed of one spinning black 
hole and one non-spinning black hole of any mass ratio. 



The complicated three-dimensional motion can be 
beautifully decomposed into two-dimensional motion in 
an orbital plane with a precession of the orbital plane 
superposed Through this modular decomposition, 
we are able to define a complete taxonomy of all three- 
dimensional orbits in terms of orbits that are closed in the 
orbital plane. Further, we find the spectrum of orbits for 
a given black hole system. Importantly, the spectrum in 
the strong-field regime shows zoom-whirl behavior dur- 
ing which an orbit sweeps out to apastron and back in a 
zoom followed by a multiplicity of nearly circular whirls 
around the center of mass. It must be emphasized that 
our results prove zoom-whirl behavior is ubiquitous even 
for comparable mass binaries and unambiguously quan- 
tified in our taxonomy. Importantly, zoom-whirl motion 
has already been observed in numerical relativity fv^ . 

From the outset, we acknowledge that the PN approx- 
imation is pushed to its breaking point in the strong- 
field regime where zoom-whirl behavior is most preva- 
lent. However, the method we advocate - locating a pe- 
riodic skeleton in an orbital basis - can be applied to 
any description of black hole binaries, including the ef- 
fective one-body (EOB) approach [l^ and extreme mass 
ratio inspirals (EMRIs) modeled by the Kerr spacetime 
[10] . Additionally, the closed-orbit taxonomy offers a new 
terrain for the comparison of the PN expansion to fully 
relativistic treatments [20|, Hi], |22[ . Although quantita- 
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FIG. 1: Top: The orbital plane precesses around the J = k 
axis through the angle ^. Bottom: The orbital plane can be 
spanned by the vectors (X, Y) or the vectors (n, §). 



tive results will change in improved approximations, the 
qualitative features should be robust, as the detection of 
zoom- whirl behavior in fully relativistic numerical exper- 
iments implies [13] • 



I. PREVIEW 

To provide the reader with a road map through inter- 
mediate results accumulated on our way to the periodic 
taxonomy, we preview some highlights here. 

Our method can be broken into two main steps: 
1. Simplified Equations of Motion. Since the 
perihelion precesses and the orbital plane precesses, the 
motion around a spinning black hole depends on angles 
as well as on the radius. In usual spherical coordinates, 
the equations of motion are quite complicated. By work- 
ing in a non-orthogonal orbital basis, we show explicitly 
that the equations of motion are independent of (non- 



orthogonal) angular variables. 

Physically, we exploit the observation^ that the orbit 
lies in the plane spanned by the coordinate r and its 
canonical momentum p; that is, the orbital plane is per- 
pendicular to the orbital angular momentum, L = r x p. 
The plane itself then precesses around the constant to- 
tal angular momentum, J. The importance of the or- 
bital plane was clear in some of the earliest papers on 
spin-precession [3], although that early work generally 
imposed a quasi-circular restriction on the orbits. We 
decompose all motion into precession of the perihelion 
within the orbital plane with a precession of the entire 
plane superimposed, with no restrictions or approxima- 
tions. A preview of the explicit construction is shown 
in Fig. [TJ A fully precessing orbit is shown on the top 
in Fig. [2] while on the bottom the orbital plane traps a 
much simplified orbit, reminiscent of the equatorial orbits 
of Kerr black holes [1^ . 

The simplified Hamilton's equations immediately in- 
form us that all eccentric orbits have constant aphelia 
and perihelia.^ When the aphelia and perihelia are one 
and the same, we have non-equatorial constant radius or- 
bits, also known as spherical orbits (as previously found 
in [H, [2^ ) . The spherical orbits are not necessarily peri- 
odic; they fill out a band on the surface of a sphere. They 
are nonetheless significant in our campaign to fully dis- 
sect the dynamics and are treated in a companion paper 

The simplified Hamilton's equations show that zoom- 
whirl patterns will be symmetric from one radial cycle 
to another when viewed in the orbital plane, as can be 
seen on the right of Fig.[2l A related subtle feature is that 
the three coordinate velocities in the orbital basis depend 
only on radius and are therefore periodic as an orbit ex- 
ecutes a radial cycle from apastron to apastron. Taken 
together these symmetries in the orbital plane are in- 
triguing for gravitational wave analysis. The waveforms 
must be decomposable into the orbital basis and there- 
fore Fourier decomposable into the three fundamental fre- 
quencies that are the time average over one radial cycle 
of the instantaneous velocities. It remains to be seen 
how advantageous this might be for gravitational wave 
astronomy. 

We restrict ourselves to completing the dynamical pic- 
ture in this article since it is the dynamics that shapes 
the gravitational waves. In this spirit, step 1 above allows 
us to proceed to step 2: 

2. Taxonomy of Fully 3D Orbits. We offer a 
method to completely taxonomize the dynamics with the 
restriction that only one of the black holes spins. Our 
approach includes all fully three-dimensional orbits de- 



^ The orbital plane is also emphasized in applications of PN dy- 
namics to pulsar timing p3l . [Sj. [25I . [26| . 

^ The constancy of periastron and apastron for every orbit might 
must have been implicitly understood in Refs, 
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FIG. 2: Top: Fully three-dimensional orbit. Bottom: The 
trajectory as captured in the orbital plane. 

scribed by the third-order Post-Newtonian (3PN) Hamil- 
tonian plus spin-orbit couplings. 

Our taxonomy extends the periodic tables for Kerr 
equatorial orbits [I^ to fully non-equatorial orbits of 
comparable mass black hole binaries. In Ref. [l9|], we 
introduced a taxonomy for equatorial Kerr motion with 
the following salient features. Each entry in the Kerr pe- 
riodic tables of Ref. [19|] is a perfectly closed equatorial 
orbit identified by a rational number 

q = w^- (1) 

z 

where w counts the number of whirls, z counts the num- 
ber of leaves, and v indicates the order in which the leaves 
are traced out. Since the rationals are dense on the num- 
ber line, the periodics are dense in phase space. Con- 
sequently, any generic equatorial orbit can be arbitrarily 
well- approximated by a nearby periodic orbit. In this 
way, any generic orbit is approximately equivalent to a 
high- leaf orbit (high z). Additionally, any generic orbit 
can be approximated as a precession around a low-leaf 
orbit, a technique that might ultimately benefit signal 
extraction. 

Our ambition in this paper is both to extend the tax- 
onomy to comparable mass binaries and to resolve the 
non-equatorial motion of spinning binaries. Truly peri- 
odic three-dimensional motion follows when the trajec- 
tory closes in the orbital plane and the precession of the 



entire plane closes simultaneously. Fully closed motion 
requires two rationals, each representing a ratio of fun- 
damental frequencies. And although in principle there 
must exist orbits that are fully periodic in the three- 
dimensional motion - as Poincare argued [s^ - our tax- 
onomy of bound orbits needs only the weaker condition 
of periodicity in the orbital plane. Not every orbit that is 
closed in the orbital plane will be closed in the full three- 
dimensional space. In other words, for the less restrictive 
condition of orbital plane periodicity we only need one ra- 
tional ratio of frequencies. As will be explained in detail 
in gilll the aperiodic orbit of Fig. [2] can be approximated 
as a precession around a 4- leaf clover in the orbital plane. 

Although the PN approximation is poor in the strong- 
field, the qualitative results should survive a full relativis- 
tic treatment. Spin-spin couplings will impose additional 
modulations on the orbital plane picture but since spin- 
spin couplings are higher-order in the PN expansion, the 
expectation has been that their effect can be treated as 
a perturbation [28] . In a companion paper we will argue 
that although a small perturbation, the spin-spin cou- 
plings are responsible for the emergence of chaos around 
unstable orbits 0. 

We emphasize that for real spinning astrophysical 
black holes, the orbits we resolve in this paper are not an 
exotic subset of orbits, but rather are descriptive of all 
bound orbits - all non-circular orbits are captured in the 
spectrum of rationals. There is a long-standing argument 
that black hole binaries will circularize by the time they 
enter the bandwidth of the gravitational wave observa- 
tories. However this is not possible for spinning black 
holes. Circular orbits do not exist for misaligned spins. 
Although spherical orbits do exist, they are destroyed by 
the spin-spin effects [29|. What's more, black hole pairs 
formed in dense clusters are not expected to circularize 
by the time of merger and are expected to be plentiful 
sources for advanced LIGO [3l[ While we restrict our- 
selves to one spinning black hole and one non-spinning 
in this paper, the scenario is both astrophysically pos- 
sible in its own right and theoretically important to lay 
the foundation for the two spinning case with spin-spin 
included, a task we return to in a companion paper [2y] 

We express Hamilton's equations in a non-orthogonal 
orbital basis in ^ We discuss the closed orbit taxonomy 
in ^IIII In ^IV( we show periodic tables for two different 
black hole binaries, a comparable mass binary and a non- 
spinning extreme mass ratio pair. Appendix €K\ details 
the projection of Hamilton's equations onto our orbital 
basis. In the conclusions, ^Vj we discuss the modulations 
predicted from spin-spin couplings and those imposed by 
spinning both black holes. 



II. HAMILTON'S EQUATIONS OF MOTION IN 
THE ORBITAL BASIS 

The culmination of this section will be the compact 
form of the equations of motion (Eqs. (|28|) ) in a non- 
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orthonormal orbital basis. To get there will require a few 
short subsections. We begin with the 3PN Hamiltonian 
including spin-orbit couplings. 

A. The 3PN Hamiltonian + Spin-Orbit Couplings 

Although the 3PN Hamiltonian is nearly a page long, 
the Hamiltonian formulation of black hole pairs has cer- 
tain advantages over the Kerr fully relativistic descrip- 
tion of test particle motion around a single black hole. 
Most notable in this context, the ADM-Hamiltonian ef- 
fectively describes center of mass motion in flat space. 
This will allow us to manipulate spatial vectors at will 
and locate the orbital plane. (This work has suggested a 
means to generalize to the fully relativistic Kerr system. 



where 



The reduced spin-orbit Hamiltonian is 

Hso = S^'^ (6) 
where in this paper, we restrict our analysis to only one 



That research is in progress.) 

To begin, take the 3PN Hamiltonian including spin- 
orbit coupling as it is conventionally written in dimen- 
sionless coordinates. If TZ is the ADM coordinate vector 
and V is the ADM momentum vector, then the dimen- 
sionless center-of-mass coordinate vector is r = and 
its canonical momentum isp = V/ /a where the total mass 
is M = mi +m2 for a pair with black hole masses mi and 
m2 and the reduced mass is /i = 17111712/ M with the di- 
mensionless combination 77 = ji/M. All vector quantities 
will always be in bold so that r is to be understood as 
the magnitude r = v^iw*. Unit vectors such as n = r/r 
will additionally carry a hat as well as being bold. The 
dimensionless reduced Hamiltonian H = H/ 11^ where H 
is the physical Hamiltonian, can be written to 3PN order 
as [32|,[33|,[3i,[3i,[3i,[33 



(2) 



(3) 



(4) 



(5) 



spinning body with reduced spin 



S = a(m?//iM) , (7) 



H = HpN + Hso , 



HpN = Hn + HiPN + H2PN + HspN 



HiPN 



{in -I) (p^)'-^[(3 + r,)p^+,?(n.p)^] ' 



r 2r^ 
2 



H2PN = ^ (1 - 57, + bif) (p')' + g [(5 - 207, - 37?2) (p2 
-27?2(fi.p)V-37?2(fi.p)4]i 
+ ^ [(5 + 87,) p2 + 37,(n . p)2] 1 _ 1 (1 + 37?) 1 



H^PN = (-5 + 357, - IQif + i^if) (p= 



.4 1 



(-7 + 42/^ - b?>rf - brf) (p^ 

1 

r 



+ (2 - 37,)7,2(n . p)2(p2)2 + 3(1 _ ^)^2(f^ . p)4p2 _ ^^Z^^ . p)6] 

^{-21 + 1367, + 1097,2)(p2)2 + 1(17 + 307?)7?(n • p) V + :i^(5 + 437,)7,(fi • p)^ 
16 16 12 

{ ik + ^^^^ ~ ^^"^^^ ~ ^^^'^'^ " 6i ^^"^^ + + ^^^''^ "^^^ ' p^'} 

^[l2 +(872 - 637r2)7,]l . 
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mass mi, and 



^1 



3m2\ 
2rrnJ 



V 



(8) 



Physical values of the dimensionless spin amplitude range 
over < a < 1. In a conipanion paper, we will generalize 
to two spinning bodies [2^. We omit spin-spin coupling 
terms. The reduced orbital angular momentum is 



r X p 



(9) 



Notice that with units included, the physical orbital an- 
gular momentum is L/iM. 

The equations of motion are given by 



dH 
dp 



dH 
dr 



(10) 



and the evolution equation for the spins and the angular 
momentum can be found from the Poisson brackets: 



B ]-f 

S = {S,H} = — xS 



(11) 



which comes to 



S = .5i 



L X S 
S X L 



(12) 



B. Conserved Quantities 



It is well-known that this system has many useful con- 
served quantities.^ The Hamiltonian is conserved by con- 
struction. The conservation of total angular momentum 
follows from Eqs. ([12]) 



J = L + S 



(13) 



Also conserved are the magnitude of S and L, as can be 
confirmed by taking the dot-products with Eqs. (fT^ : 

S.S=i^(S'2)(xL-(SxL)=0 

L.L = i|(L2)«S-(LxS)=0 (14) 



^ Although it is by now well-confirmed that there is chaos when 
the black holes spin p^, [H, HO, |4l|, IH, El, |3 , we are deal- 
ing with a restricted situation of only one body spinning in the 
Hamiltonian formulation and our orbits are not chaotic to this 
order in the approximation [45[ . This is not in confiict with ear- 
lier work on chaos in the Lagrangian approximation [46[ . At 
higher order including spin-spin couplings, a pair of spinning 
black holes loses constants of the motion opening a window for 
chaotic motion 1291. 



Finally, the component of L in the J direction must be 
conserved as can be seen from 



L = Si 



S X L 



J X L 



(15) 



from which it follows that the change in L is always per- 
pendicular to J. 



C. The equations of motion 

We want to express the equations of motion derived 
from Eq. (p!Q|) in the following form: 



r 
P 



Ap 
Cp 



Bn- 
Dh 



spm pieces 
- spin pieces 



(16) 



This form helps consolidate the equations of motion be- 
fore we project from the vector equations to component 
equations in the next section ( ^IIPp . We need to iden- 
tify the functions A^B^C^D in terms of derivatives on 
the Hamiltonian, which can be thought of as a function 
of (r, p, (n-p)). Considering the non-spinning piece, Hpn 
first, we break up the partial derivatives on the right hand 
side of Hamilton's equations in the following way. 



dH 



PN 



dp 



dH^ 



PN 



0>p2 



dp^ 

r,(fi-p) 



dH 



PN 



5(n • p) 



g(n • p) 
dp 



and 



dH 



PN 



dr 



dH 



PN 



d{h-p) 



d{n • p) dh dHpN 



r,p 



dh dr 



dr 



P,(np^ 



dr 
d^' 



where we are careful to indicate the quantities held fixed 
in each term. Using 



d{n • p) 9n p (n • p) ^ 

dh dr r r 



we define 



A = 2 
B 



dH^ 



PN 



9p^ 

dHpN 



(17) 



(18) 



•,(n-p) 



r,p 



c 

D 



d{h-p) 

1 dHpN 
r d{h • p) 



dH 



PN 



dr 
dHpN 



B 

r 

dHpN 



dr 



P,(n p) 



P,(n p) 



(n-p) 



d{n ■ p) 



(ft-p)C . 



(19) 
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The variations of the spinning piece of the Hamiltonian 
are simply 



dHso 

dp 
dHso 



Si 



S X r 



-Si 



p X S 



, L-S 

Ml — 3— n 



(20) 



dr r"* 

With these definitions, and making use of— pxS = Sxp, 
we can write the vector equations of motion compactly 
as 

S X r 



r = Ap + 5n 
p = Cp^ Dh 



Sxp 



3Si 



L S 



-n 



(21) 



To go from these vector equations to component form 
requires we choose a basis. In the next section we will 
build the orbital basis of Fig. [H and cast Eqs. ([2T]) in 
component form. 



where Ol = arccos(L • J) is constant. The X axis is or- 
thogonal to both J and L by construction and so lies on 
the intersection of the orbital plane and the equatorial 
plane. The entire orbital plane maintains a fixed angle 
Oy = 7t/2 — Ol with J as it precesses with the precession 
of L. The motion of the orbital plane can be understood 
as the motion of X in the equatorial plane through an 
angle ^ where 



^ = J X X 



(24) 



Let J = k and let the unit vectors i and j span the 
equatorial plane. The speed of ^ motion can then be 
determined: 



cos ^ = X • i = ^.^l'^ = Ql^ ■ i = -Ql sin ^, (25) 



D. The Orbital Basis 



where we have explicitly used the constancy of the mag- 
nitude L. From this we conclude 



The clarity of the form of the equations of motion de- 
pends on the basis used to express them. There are sev- 
eral choices although the one we are calling the orbital 
basis leads to profound clarity of expression. We will 
build a non-orthogonal, unit normalized basis (n, ^, ^) 
in this section. 

There are two special planes to consider when spin 
precession drives three-dimensional orbits. There is the 
orbital plane, which is the plane perpendicular to L, and 
there is the equatorial plane, which is the plane perpen- 
dicular to J (see Fig. [1]). We will find an orthonormal 
basis that spans the orbital plane and then add the mo- 
tion of the plane itself, which will be in a direction that 
is not orthogonal to the orbital plane. The technique of 
moving into an orbital plane and projecting equations of 
motion onto this basis is familiar from celestial mechan- 
ics and has see n appli cation in the PN approximation to 
binary pulsars [U |25|, [2^ . We depart from the usual 
approach by adopting a non-orthogonal basis. 

The vectors n and p lie in the orbital plane by the 
definition L = rn x p. We can also span the orbital 
plane by orthonormal vectors (fi, ^) where 



L X n 



(22) 



We will work in terms of (fi, ^) when considering motion 
in the orbital plane. 

To separate out the precession of the orbital plane from 
the three-dimensional motion, another basis will be use- 
ful for intermediate steps. The (X,Y) basis spans the 
orbital plane but rotates with the precession of L: 



X 
Y 



J X L 

sin Ol 
L X X 



(23) 



nL = ^ = Si- 



(26) 



does not depend on any angles. 

To find our simplified equations of motion we work in 
the non-orthogonal, unit normalized basis (fi, ^, ^) in 
the next section. 



E. Final equations of motion in the orbital plane 

The four equations of motion in the orbital plane are 
obtained by projecting Hamilton's equations onto the ba- 
sis vectors (fi, ^). We do this explicitly in appendix ^Al 
where the projections onto the orbital basis vectors gen- 
erate the four equations. 



dH 



r # 



dp 

dH . 

• ^ 

dp 

dH 



pn= -— - n 
or 

p.$ = — — ^ 
or 



(27) 



Compiling the equations of appendix ^ concisely gives 
the remarkably simple equations of motion in the or- 
bital plane coordinates (r, <l>) and and their canonical 
momenta (P^,P<|,), 
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n ST 


r = APr 


+ B, 


Pr = - -P^ + D + 3,5i— ^ 














P<j, = 



(28) 



where = L. The orbital plane precesses at a rate that 
depends only on r, as calculated in the previous section 

{mm- 

^ = nL = Si^ , P^ = (29) 

with Si = (^2 -\- T] and = Lz. For completeness, 

we can track the precessions of the spin and the angular 
momentum: 

S I^lJ X S 

L = I^lJ X L . (30) 

As noted earlier, S • L =constant. Notice, this basis is 
explicitly constructed for S x L 7^ 0. When the spin and 
orbital angular momentum are aligned, anti-aligned, or 
spin is zero, then motion is confined to a plane and we 
should use the usual equatorial planar basis. This is done 
explicitly in ^TVBl 

The functions A, 5, C, D of Eqs. (p!9|) depend only on 
r, Pr and constants. This can be seen by noting that 
^4, 5, C, D are functions of (r, (fi • p), p^). Now, p can be 
written in terms of a piece in the radial direction and a 
piece perpendicular to the radial direction: 

p = (p • n)n + (n X p) X n 

= P^n+-xn (31) 
r 

so that 

P' = Pr+^ (32) 

and L equals a constant. The term p^ can therefore 
be expressed as a function of (r, Pr) only. Meanwhile, 
(n • p) = Pr and so any function of (r, (fi • p),p^) is 
equally well a function only of (r, Pr) and constants. 

Consequently, the above equations of motion are, 
amazingly enough, independent^ of the angles (^,^). 



^ The simplicity of the orbital equations of motion can of course 
be recast in terms of symmetries. Rotations about L through 
the angle ^, and rotations about J through the angle 

^, leave the dynamics invariant. The coordinates $ and 

^ are cyclic and = P^^ = 0. The symmetries correspond to 
conserved L and Lz with = L and = Lz. 



The purely radial dependence of the equations of mo- 
tion immediately informs us of four crucial facts (valid 
to this order in the PN approximation): 

1. There are constant radius orbits. These were al- 
ready found in Ref. [H, [2^ and drop out particularly 
simply in the orbital basis. In the equatorial plane these 
are of course the usual circular obits. Out of the equa- 
torial plane they have been called spherical orbits in the 
literature since the orbits trace out an annulus on the sur- 
face of a sphere. We will continue in this spirit and call 
them spherical orbits, or, more exactly, constant radius 
orbits. To see that they exist, notice that the spherical 
obits correspond to solutions of r = = and these 
roots - according to Eqs. ([28|) - can only depend on con- 
stants (mi/m2, S^L^Ols)^ not on angles. In fact, we have 
a pseudo effective-potential description as we explain in 
the next section. 

2. All eccentric orbits have constant apastra and per- 
ihelia. Similar to the reasoning above, the solutions to 
the condition = are the turning points^ and thus are 
the apastron and perihelion of any orbit. These values 
can depend only on constants [2^. 

3. The orbital angle and precessional angle swept 
out between successive apastra are constant. The angles 
swept out as an orbit moves from one apastron to another 
are simply given by 

A$ = 2 / -dr 

= 2 / -dr (33) 

and depend only on constants. 

4. All three coordinates velocities are periodic in r. 
The three coordinate velocities in the orbital basis 

r , <i> , and ^ , (34) 

depend only on the variable r (or Prij)) and thereby 
inherit r's periodicity. They can be averaged over one 
radial cycle to define three fundamental frequencies 

uOr = — , ^cE) , and cj^ (35) 



^ When P^ = 0, P = and so f = APr + P = 0. Therefore we 
take the f = condition to be synonymous with Pp = 0. 
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FIG. 3: A generic orbit in the strong- field. The initial con- 
ditions are 7712 /mi = |, L = 3.5, Ols = f , a = 0.9, and 
Ti = 18. Top: The full three-dimensional orbit. Middle: A 
projection of the full orbit onto the equatorial plane. Looking 
closely, the angle swept out from leaf to leaf is not the same 
under this projection Bottom: The orbit as caught by the 
orbital plane. The angle swept out in the orbital plane from 
leaf to leaf is always the same. Further, the constancy of the 
apastron is clear. 



where Tr is the radial period from apastron to apastron 
(We will define uj^ and cj^ a bit later in Eqs. ([36|) and 

(EID).^ 

The powerful simplicity of the description in the orbital 
basis is visually manifest in Fig. [3l The figure shows a 
fully three-dimensional, generic orbit on the top. The 
lowest panel is the same snapshot captured in the or- 
bital plane. Notice how the orbital plane reveals the con- 



Points 1 and 2 are also true for Kerr non-equatorial orbits, al- 
though 3 and 4 are not. 
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FIG. 4: The initial conditions are 1712 /mi = |, L = 3.5, 
Ols = ^,a = 0.9 and E = -0.023548373360051289666. Top: 
The full orbit. Middle: A projection on the equatorial plane. 
Bottom: The orbit as it appears in the orbital plane 



stancy of the three-dimensional apastron and perihelion 
as claimed in point 2 above. Also notice that the spacing 
between leaves is always symmetric in the orbital plane. 
Said another way, the angle swept out in the orbital plane 
between apastra is always the same, as claimed in point 
3. Neither of these features is apparent from the fully 
three-dimensional snapshot or from the projection onto 
the equatorial plane shown in the middle view. Another 
generic orbit is shown in Fig. [H 

The four facts above have two significant implications: 

• A gravitational waveform can be Fourier decomposed 
in the three fundamental frequencies of Eq. . 

• There should be a spectrum of orbits that are closed 
in the orbital plane, and that spectrum must have a cor- 
respondence with the rationals. We can therefore gen- 
eralize the Kerr taxonomy of Ref. [l^ to non-equatorial 
orbits of comparable mass black hole binaries. 
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We complete the dynamical picture by moving on to 
the periodic taxonomy for black hole binaries. 



III. CLOSED ORBIT TAXONOMY 

Any dynamical study benefits from locating the closed 
orbits - orbits that return to their initial values after 
some finite period. Poincare was the first to reajize that 
periodic orbits structure the entire dynamics [30|. Al- 
though a set of measure zero, the periodic orbits are dense 
in phase space. Consequently, any orbit can be approx- 
imated as near some periodic orbit. In that sense, the 
periodic set forms the skeleton of the dynamics. What's 
more, they are all one needs to know since to arbitrary 
precision even an aperiodic generic orbit is arbitrarily 
close to some periodic orbit, though possibly one with 
very high period. 

The periodic set corresponds to a spectrum of ratio- 
nal numbers. That spectrum of rationals shows that the 
zoom- whirl behavior known for extreme-mass-ratio bina- 
ries is prevalent in the strong-field regime of comparable 
binaries as well. Zoom- whirl behavior is therefore not ex- 
otic but rather the norm for the strong-field. The spec- 
trum of rationals renders the zoom- whirl behavior of any 
orbit quantifiable and unambiguous. 

Consider the coordinate velocities of Eq. ([34|) . Taking 
the time average of the <l>-frequency over one radial cycle 
gives the fundamental frequency 



so that 




An orbit that is closed in the orbital plane has rationally 
related frequencies 

^ = 1+9* . (38) 

UJr 

By Eq. ([37j) . we can interpret the rational in terms of the 
angle swept out from leaf to leaf in the orbital plane 

— = l + g$ = l + ^<i, + — , (39) 

where we have written the rational in terms of a triplet 
of integers, as can always be done [19]. In the equatorial 
case we know that A<I> > 27r for all eccentric orbits; all 
relativistic orbits overshoot as the famous precession of 
the perihelion of Mercury attests. For that reason we 
have separated out a 1 from the definition of in Eq. 
(|39l).^ 



^ There is an interesting anomaly that has to be mentioned. As it 
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FIG. 5: A = 1/3 orbit. The initial conditions are 7x12/ mi — 
\,L = 3.5, Ols = f , a = 0.9, and E = -0.0220582156. Top: 
Fully three-dimensional orbit. Bottom: The trajectory is a 
3-leaf periodic in the orbital plane. The first radial cycle is in 
bold. 



By analogy with the equatorial Kerr case of Ref. [l^ , 
counts the number of leaves (or zooms) , specifies 
the order in which the leaves are traced out, and 
counts the number of additional full 27r whirls taken be- 
tween apastron and apastron. To clarify the role of v^^ 
label the leaves sequentially through z — 1 starting with 
the initial apastra. Then equals the number of the leaf 
that the orbit jumps to after the starting apastron. The 
meaning of the rational is best illustrated with an exam- 
ple. An orbit with g<i> = 1/3 is shown in Fig. [51 This is 
a 3-leaf orbit {z^ = 3) that moves to the first leaf in the 
pattern {v^ = !)• Since = 0, there are no additional 
whirls from leaf to leaf. In Fig. [6l we show a g.^, = 2/3 
orbit. That is, a 3-leaf orbit {z^ = 3) that moves to the 
second leaf in the pattern {y^ = 2). Since = there 



happens, in this 3PN approximation, it is possible for periodics 
in the orbital plane to undershoot 27r; that is, from one apas- 
tron to another A<l> < 2tt < 0). This never happens with 
Kerr equatorial orbits and may just be a peculiarity in the ap- 
proximation, although we won't know for certain until the Kerr 
non-equatorial case is completed. Also notice that while some 
orbits may undershoot in the orbital plane, $ is not the whole 
story (see the footnote in ^IV B|) and some of the apparent re- 
gression is more than compensated for by ^ . 
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FIG. 6: A = 2/3 orbit. The initial conditions are 1712/ mi = 
^, L = 3.5, Ols = f , = 0.9, and E = -0.0211669686. Top: 
Fully three-dimensional orbit. Bottom: The trajectory is a 3- 
leaf periodic in the orbital plane that skips a leaf each radial 
cycle. The first radial cycle is in bold. 
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FIG. 7: An orbit for which = The initial condi- 

tions are 7112/ mi = |, L = 3.5, Ols = f , a = 0.9, and 
E = -0.0220323426 Top: The full three-dimensional orbit. 
Bottom: The orbit in the orbital plane is a precession of the 
exact = 1/3 orbit. 



are no additional whirls from leaf to leaf. For a z^-ledif 
orbit, the range of for orbits that overshoot, that is 
precess, is 



1 < < — 1 ^ 



if > 1 
if z^ = 1 



(40) 



To avoid degeneracy, we require that z^ and be rel- 
atively prime, or in other words, that = 2/4 is the 
same diS q^ = 1/2. 

The periodic orbits in the orbital plane are a set of 
measure zero in the space of orbits, just as the rationals 
are a set of measure zero in the set of the real numbers. 
However, just as the rationals are dense on M, the peri- 
odics are dense in the space of orbits and so any generic 
orbit can be arbitrarily well approximated by an orbit 
that is closed in the orbital plane. 

For instance, the orbit of Fig. [3 is very near the 3-leaf 
orbit and can be interpreted as a precessing 3-leaf orbit. 
Or we could do better by approximating this orbit as a 
q^ = 67/200, that is, an orbit with 200 leaves that skips 
to the 67th successive leaf in the pattern with each radial 
cycle. 

By the same token, the randomly selected orbit of Fig. 
[3] is very nearly di q^ = 1/5, that is, a 5-leaf clover, and 
that of Fig.[3]is very nearly a g<i> = 7/25 - an orbit with 25 



leaves that skips to the 7th successive leaf in the pattern 
with each radial cycle. 

It is important to notice that although the orbit of 
Fig. [5] closes in the orbital plane after 3 radial cycles, it 
does not close in the full three-dimensional space since 
the orbital plane has not returned to its original location 
after only 3 radial cycles. A fully closed orbit also has to 
close in ^. Taking the time average of the rate of change 

over one radial cycle gives the fundamental frequency 



so that 



cj^ _ A^ 

COr 27T 



(42) 



The average precessional frequency may not be rationally 
related to the radial frequency for a rational q^: 



(43) 



where by we mean any real number, not just a ra- 
tional. This time we do not separate out a 1 from the 
definition of the number. So, represents the fraction 
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of 27r swept out as the plane precesses. The orbit of Fig. 
[5] has a cr^ ~ 0.346..., where we have only listed the first 
3 significant figures. Although numerical imprecision of 
the computer truncates this at a finite number of digits, 
and therefore effectively approximates by a rational, 
it is in principle an irrational. After 3 radial cycles, the 
orbit has closed in the orbital plane but not in three di- 
mensions. The entire orbital plane has overshot its initial 
location by 3cr^ — 1 0.038.... Therefore, even the 3-leaf 
clover in the orbital plane of Fig. [5] will fill out the surface 
of the three-dimensional picture. 

In fact, from the equations of motion Eqs. ([28|) - ([29|) . 
we know that 

l^q^ = -jCT^^f{E,L) , (44) 

where f{E, L) is a function of L through the depen- 
dence of the apastra and periastra's dependence on E^ L. 
This curve for (m2/mi = 1/4, a = 0.9,^^5 = tt/S^L = 
3.5) is shown on the top of Fig. [HI Eq. (|44|) is interest- 
ing for two reasons. It means that is not generally a 
rational number when is rational. What's even more 
interesting is that this does not seem to matter. Because 
effectively fixes the value of for a given L through 
relation (j33j) , if a generic orbit is well approximated by a 
periodic in the orbital plane, the precession of its orbital 
plane is also nearby in the phase space sense. ^ 

This point is emphasized on the bottom of Fig. [HI 
which shows the circle traced out by the precession of 
the L vector for the 3-leaf clover in the orbital plane of 
Fig. [5l The straight line from the origin to the ring indi- 
cates the direction of L after 3 radial cycles have elapsed 
and the 3-leaf clover has executed one complete period 
in the orbital plane. The of this orbit is ~ 0.346... so 
that in 3 radial cycles the orbital plane has just overshot 
its initial location. By comparison, the nearby orbit of 
Fig. [7] has di ^ 0.348... and its orbital plane similarly 
has just barely overshot its initial location. The orbit of 
Fig. [7] precesses around the 3-leaf clover of Fig. [5l and 
its entire orbital plane precesses around J, sticking close 
to the precession of the 3-leaf clover's orbital plane. In 
fact, the precessions of L are superposed on the bottom 
of Fig. [8] and the difference between them is impercepti- 
ble at the resolution shown. The fact that the precession 
of L for these two orbits is effectively indistinguishable 
confirms that the two orbits are not only near each other 
in the orbital plane, they are genuinely near each other 
in 3D as well. 

Formally, the above argument ensures that any orbit 
can be approximated as arbitrarily near an orbit that 



^ Of course, we can alway approximate any irrational, including 
cr^ by a rational. But then we are describing approximately 
periodic orbits as opposed to orbits that are formally exactly 
periodic and there doesn't appear to be any advantage in taking 
this tack. Hereafter, we will consider cr^ to be generally irrational 
for any rational q^. 




"^"^'"^^---^^^ 





FIG. 8: (m2/mi = 1/4, a = 0.9, Ols = 7r/3,L = 3.5). Top: 
versus cr^. The dots mark (q^ = 1/3, cr^jf ~ 0.346...) and 
{q^ — 2/3, cr^ ~ 0.679...). Bottom: The circle traced out by 
the tip of the vector L for the orbit of Fig. [5] The straight line 
represents the L vector when 3 radial cycles have elapsed and 
the q^ = 1/3 orbit has closed in the orbital plane. The same 
plot for the precessing 3-leaf clover of Fig. [3 is superposed al- 
though the two are so close that they cannot be distinguished 
in the graph. The fact that they cannot be distinguished con- 
firms that the two orbits are genuinely near each other in 3d 
as well as in the orbital plane. 



is periodic in the orbital plane with the same L. In 
other words, the orbital periodic spectrum for a given 
L is dense. If we remove the restriction of comparing 
orbits of the same L, it follows that the set of orbits pe- 
riodic in the orbital plane is dense in the entire space. 
The argument can be sketched as follows. According to 
Poincare, the set of orbits that is fully periodic in 3D is 
dense in the phase space. This set is a subset of the or- 
bital plane periodic set. Therefore, if the subset is dense, 
the set itself must be dense. 

In short, we can understand the entire three- 
dimensional orbital dynamics through orbits that are 
closed in the orbital plane and that one rational, not 
two, is needed for a taxonomy. These conclusions are of 
course only valid up to 3PN with spin-orbit coupling. In 
the summary we will discuss the modulation expected by 
going to higher order in the approximation. In the mean- 
time, we move on to the periodic tables for comparable 
mass binaries. 
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FIG. 9: (m2/mi = 1/4, a = 0.5, Ols = 7r/4). The pseudo 
effective-potential as a function of r shows a maximum at 
the unstable spherical orbit and a minimum at the stable 
spherical orbit. The homoclinic orbit is indicated with a 
dashed line. The other lines correspond to, in ascending or- 
der, = f , -I, §• The higher orbits quickly stack together 
near the homoclinic orbit. 



IV. PERIODIC TABLES 

Since every orbit can be approximated as one that is 
periodic in the orbital plane - just as every irrational 
can be approximated by a rational - we can build a ta- 
ble of orbits for black hole binaries with a given mass 
ratio, spin, and angle between the spin and the orbital 
angular momentum {m2/rrii^ L^Ols)- Such a periodic 
table works in analogy to the chemical periodic table as 
illustrated in Fig. [121 Every entry represents a closed 
orbit labeled by a rational. The energy and the rationals 
both increase monotonically from top to bottom and then 
from left to right. Unlike the chemical periodic table, the 
black hole periodic tables are infinite since the rationals 
are infinite and we show only a handful of entries. 

When we discuss specific tables, we will always 
take all entries in a given table to have the same 
(m2/mi, 5*, L^Ols)' The entries vary only in energy (and 
therefore in q^). Periodic tables can be constructed for 
any values of these parameters. However, some ranges 
give fuller tables in the sense that certain ranges permit 
more whirliness. To understand the ranges of whirliness 
requires careful consideration of properties of the spher- 
ical orbits. Several results concerning non-equatorial 
spherical orbits are presented in [2^. Here we summa- 
rize the pertinent results for constructing periodic tables 
with use of a pseudo effective-potential picture. 

A sensible condition for an effective potential formula- 
tion is that the Hamiltonian depend only on r and con- 
stants. Generally, the Hamiltonian of Eqs. ([2])- ([5]) does 
not admit a simple effective potential formulation since it 
is a complicated function of . We have already argued 
that (r, p, S) can be written as a radial Hamiltonian 
i^(r, Pr), yet it remains a complicated function or P^. 
However, if we only consider 

Kff = i^(Pr = 0) , (45) 



then we have a good representation of the effective po- 
tential at the turning points. We cannot misuse the V^ff 
by trying to interpret motion away from the turning 
points, but it gives a perfectly valid description of the 
behavior at aphelia and periastra as well as on spheri- 
cal orbits. Hereafter we'll shorthand the term "pseudo 
effective-potential" by "effective potential" . An effective 
potential for comparable mass binaries is shown in Fig. 

Evident is the lowest energy orbit at the stable spher- 
ical orbit. The highest energy non-plunging orbit is the 
unstable spherical orbit. We are interested in energet- 
ically bound orbits here, i.e. orbits with E < {). If 
the unstable spherical orbit has energy £^ < 0, then 
the spectrum of periodic orbits densely fills the energy 
range between the stable and unstable spherical orbit. If 
instead the unstable spherical orbit has energy £^ > 0, 
then the the spectrum of periodic orbits densely fills the 
energy range between the stable spherical orbit and the 
marginally bound orbit at £^ = 0. The energy levels of a 
few periodic orbits are indicated by solid lines in Fig. [9l 

Since q^ is monotonic with energy, the entries in the 
periodic table are bounded: 

(46) 

The value of qmin is the q^ of the stable spherical orbit 
and qmax is set by the q^ of the maximum energy bound 
orbit. 

The value of qmax depends on the largest energy orbit 
allowed for that L. When an unstable spherical orbit ex- 
ists there is always an orbit with the same E and P at a 
large radius. The maximum of Vef[ in Fig. [9] marks the 
unstable spherical orbit. Drawing a line of constant en- 
ergy across the effective potential locates the apastron of 
the orbit with the same E and P as the unstable spher- 
ical trajectory. When released from apastron, this orbit 
whirls an infinite number of times as it approaches the 
unstable spherical orbit and is formally a homoclinic or- 
bit; that is, it approaches the same invariant set in the 
infinite future and the infinite past. Although not strictly 
periodic - the homoclinic orbit never returns to apastron 
- it can be considered the infinite winding limit of the 
1-leaf periodic orbits 0. As such it is the = oo 
limit and we assign homoclinic orbits di q^ = oo. Conse- 
quently, if the range of parameters has a bound unstable 
spherical orbit, then it has a {q^ = oo) homoclinic or- 
bit and the associated periodic tables will exhibit much 
whirliness. 

Since, for a given P, the stable spherical orbits bound 
the allowed energy range from below, they also bound 
the value of g$ from below. One might presume that 
Qmin = but, importantly, this is not the case. To see 
this notice that a spherical orbit obviously does not have 
a radial cycle. The qmin set by the stable spherical orbit 
can instead be thought of as the value of q^ for a nearby 
eccentric orbit in the limit that the eccentricity vanishes: 

Qmin ^lim e^O 1 • (^'^) 
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This allows us to derive the Qmin from a stability analysis 
since the limit of zero eccentricity is effectively the limit 

of constant radius, which implies <i> is constant and 

that cOr = is given by a small perturbation around 
the stable spherical orbit. We then have, just as we found 
in Ref. [la] that^ 



^{rs) 



(48) 



We choose for the sake of illustration to consider ta- 
bles capable of probing high whirliness. For this reason 
we stay within the range set by the last stable spheri- 
cal orbit (LSSO) and the homoclinic orbit Llsso < L < 
Lhomociinic Of coursc, the drawback is that we are push- 
ing the PN expansion to the breaking point. Although 
these inner strong-field orbits probe beyond the confi- 
dence of the 3PN approximation, the general method of 
constructing periodic tables in a full relativistic treat- 
ment is robust, as proven by the Kerr demonstration of 
Ref. [H. 








A. Periodic Tables in the Orbital Plane 

The purpose of the Post-Newtonian expansion is of 
course to approximate the behavior of comparable mass 
binaries. A periodic table for a binary with m^jvfix = 1/4 
is shown in Fig. [TOl The heavier black hole has a spin 
amplitude of a = 1/2 and the angle is L • S = cos(7r/4). 
The orbits do not lie in a plane and are fully three- 
dimensional, like those of Figs. [3] and [H Each entry is an 
orbit that is periodic in the orbital plane, although not 
necessarily fully periodic. The energy and the rational of 
each entry increase from top to bottom and from left to 
right. 

Notice that the first two entries are blank before the 
appearance of the 2-leaf clover in entry 3. These are 
blank because, for this (rn2/mi^a^ S^L^Ols)^ the = {) 
and = 1/3 orbits simply do not exist since qmin is 
just above 1/3 (and Qmax = oo). We saw this before 
in the Kerr system [ly]. The implication is important. 
All eccentric orbits - for this range of parameters - show 
zoom- whirl behavior. None of them look like the slight 
precession of the perihelion of Mercury. 



^ Although most orbits precess in the orbital plane, Eq. (|48|) ac- 
tually allows for regression when , . "^f < 1. When the orbital 

I -^r I 

plane and the equatorial plane align (S x L) = 0, then regression 
seems intuitively obvious. It only means that ^ is not the whole 
story of the motion of n and we must also add in ^ to see that 
the actual orbit precesses in the equatorial plane, as we will do 
explicitly by moving to the equatorial basis in ^IV Bl We did 
however also see regression out of the equatorial plane and it is 
difficult to say whether it is a flaw in the PN approximation or 
if it will survive a full relativistic treatment. 



FIG. 10: A non-equatorial periodic table for which the heavier 
black hole spins with amplitude a = |, the mass ratio is 
1x12/ mi = J, L = 3.2 and Ols = f • All valid entries up to 
= 3 are shown. The final entry begins to show a departure 
from true periodicity as a result of numerical error. The high 
numerical precision required to keep the simulated orbit near 
a perfectly periodic one is a reflection of the tight stacking of 
high orbits near the top of the potential in Fig. [9] 

Every orbit in this system can be arbitrarily well- 
approximated by an entry in the table, because the pre- 
cessional motion of the entire plane is also effectively 
fixed by the rational q^ as the plot of Fig. [11] shows. 
If could be chosen independently of q^ for a given L, 
our conclusion would not follow. However, q^ versus 
lies on a one-dimensional curve. Once q^ is known, cr^ 
can be read off. Physically this means that an orbit that 
is very near di q^ = 1/2 will precess around the 2-leaf 
clover in the orbital plane and that the precession of the 
entire orbital plane will be very close to the precession of 
the true 2-leaf clover's orbital plane. 



B. Periodic Tables in the Equatorial Plane 

Some binary parameters will automatically restrict 
motion to the equatorial plane and these require spe- 
cial discussion. For instance, if neither black hole spins 
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FIG. 11: versus cr^ for the table in Fig. [TOl 

then the system is spherically symmetric and all orbits 
are confined to a plane: the orbital plane is the equa- 
torial plane. Similarly, if one of the black holes spins 
but the spin is aligned with the angular momentum or 
anti-aligned, then the motion will again be restricted to 
the equatorial plane. We summarize these three cases as 
S X L = scenarios. 

When S X L = 0, our orbital basis (fi, ^, ^) is not 
defined and must be replaced with the usual planar basis 
(n, (^), where (p is the usual angle measured^^ between fi 
and i. The equation of motion for cp is simply 

= <i) + ^ . (49) 

The q we must use for the equatorial plane of the non- 
spinning system is then the same as the one we used in 
Ref. [ii 

^ = l + g . (50) 

Each entry is specified by this one rational which repre- 
sents the ratio of the time averaged orbital angular fre- 
quency in the equatorial plane to the radial frequency, 
just as in the Kerr case of Ref. [l^. The rational can be 
read off the topology of the orbit as q = w -\-v/z] that is, 
the number of whirls, the number of leaves and the order 
in which the leaves are laid out fix q. 

The table of Fig. [T2l reflects a non-spinning black hole 
system with an extreme mass ratio of 1712/ mi = 10~^. 
The first 3 entries are blank since qmin is just below 1/2. 
Although we only show entries up to z«^> = 4 for < 
2, for these parameters qmax = 00. [121 Although we 
defer a detailed comparison between the 3PN spin-orbit 
system of this paper and the Kerr system, we point out 
the intriguing possibility that periodic tables could be 
used to further test the accuracy of the Post-Newtonian 
expansion [2l|. 



If we were epsilon out of the equatorial plane, we would see 
regression in the orbital plane because much of the motion is 
taken through ^. Intuitively then, those instance of regression 
that are just barely out of the equatorial plane are not that 
surprising. 




FIG. 12: A periodic table for orbits in a non-spinning system 
a = 0, extreme-mass-ratio system, ^ = 10~^. The angular 
momentum is L = 3.9 . Since a = ail orbits lie in the equa- 
torial plane. Periodic tables such as this one could be used to 
expand on comparisons with the full relativistic system. All 
valid entries up to = 4 are shown. As before, the final 
entry begins to show a departure from true periodicity as a 
result of numerical error. 



V. CONCLUSIONS 

To recap in bullet format, for comparable mass binaries 
with one spinning black hole and one non-spinning black 
hole as approximated by the 3PN Hamiltonian plus spin- 
orbit coupling, our main results are: 

1. Simplified Equations of Motion in an Orbital 
Basis 

From which we find 

• constant aphelia and perihelia for non-equatorial ec- 
centric orbits, and 

• three fundamental frequencies that depend only on 
radius. 

and 

2. Taxonomy of Fully Three-Dimensional Or- 
bits 

For which we find 
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• there exists a spectrum of closed orbits in the orbital 
plane corresponding to a subset of the rationals; 

• one rational, not two is required for an orbital plane 
taxonomy of constant angular momentum slices; 

• all orbits can be approximated as near an orbit that 
is perfectly closed in the orbital plane; and 

• zoom-whirl behavior is ubiquitous in comparable 
mass binary dynamics and entirely quantifiable through 
the spectrum of rationals. 

The first discovery we made in Ref. [l^ with our pe- 
riodic taxonomy for equatorial Kerr orbits is that pre- 
cessing elliptical orbits such as Mercury's are excluded 
in the strong- field Kerr regime - just as Keplerian el- 
lipses are excluded in any relativistic regime. Instead, 
at close separations, all equatorial Kerr eccentric or- 
bits trace out precessions of patterns best described as 
multi-leaf clovers, whirling around the central black hole 
before zooming out quasi-elliptically [13, 13 • ^^^^ 
per we have found that this same conclusion applies in 
the strong-field regime of our comparable mass black hole 
binaries, even out of the equatorial plane. Our periodic 
tables in the orbital plane show zoom-whirl behavior as 
the norm in the strong-field regime and not as the excep- 
tion. 

The further importance of the orbital dynamics lies in 
its direct imprint in the gravitational waveform [i^, . 
The waveform will necessarily reflect the features above. 
For instance, an equatorial circular orbit (neglecting 
radiation reaction) is described by essentially one fre- 
quency. By contrast, all other orbits in the strong- field 
regime generate highly modulated waveforms naturally 
described by harmonics of the 3 orbital basis frequencies, 
which in turn directly correspond to the natural frequen- 
cies of a nearby periodic orbit. 

Naturally, we should ask about the astrophysical likeli- 
ness of detecting any such orbits with either the LIGO or 
LISA observatories. Although estimates vary |5l[, stel- 
lar mass black hole pairs are currently the favored source 
for advanced ground-based dectectors and intermediate 
mass black hole pairs are considered an important ob- 
jective for space-based LISA science. It is challenging 
to definitively assess the spins and eccentricities of black 
hole/black hole binaries given the absence of observa- 
tional constraints [52]. Still, one can guess that long-lived 
stellar binaries that might collapses to a pair of bound 
black holes would circularized by the time the pair en- 
ters the strong-field due to angular momentum lost in 
the form of gravitational radiation. By contrast, for 
shorter-lived black hole binaries formed in globular clus- 
ters, the astrophysical likeliness of eccentric orbits slid- 
ing in the LIGO bandwidth is assessed to be > 30% for 
eccentricities > 0.1 in Ref. [3l[. All such binaries would 
necessarily transit near the periodic set on inspiral. Even 



if the inspiral happens too quickly to witness multiple 
executions near a low-leaf clover, the orbit can still be 
sewn together as a skip from a piece of one periodic to 
a piece of another. Finally, the spins and eccentricities 
of intermediate black hole binaries detectable by LISA 
are most difficult to predict although we should expect 
them to spend a more generous allotment of windings on 
eccentric orbits in the strong-field. 

Before closing, we have to mention the effects of spin- 
spin coupling on the orbital basis picture. The spin-spin 
correction introduces additional precessions of the spin 
and this destroys the constancy of the angle between S • L 
and generally introduces explicit angular dependence in 
the equations of motion. Interpreted as a small pertur- 
bation to the system here, the spin-spin couplings cause 
additional wobbling of the precessional motion. 

When both objects spin, the impact of spin-spin cou- 
pling can be particularly destructive. It is by now well 
documented that two spinning black holes in compara- 
ble mass binaries exhibit chaotic motion in the conser- 
vative system [13, S, EqI, EH, [13, S, SI. There are 
not enough constants of the motion to ensure regular be- 
havior. Even the restricted spin-orbit scenario of this 
paper is clearly vulnerable to chaos even perturbed as 
it admits a homoclinic orbit. Under perturbation, ho- 
moclinic orbits tend to disrupt into a homoclinic tangle, 
an infinite intersection of stable and unstable manifolds. 
At root, chaos emerges as periodic orbits proliferate in 
a bounded region of space. Our clean taxonomy of pe- 
riodic orbits corresponding to a simple set of rationals 
would give way to a glut of periodic orbits corresponding 
to some fractal set, as in systems with a strange repeller 
- the Hamiltonian analog to the strange attractor. The 
complex motion may be damped by energy and angular 
momentum loses to gravitational waves but, at the least, 
chaos signals the breakdown of the simple set of peri- 
odic orbits. The onset of chaos can be directly identified 
with the spin-spin couplings, and we leave this task to a 
forthcoming companion paper . 
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However, when spin-spin coupling is included, there are no cir- 
cular or even spherical orbits [29[. 
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APPENDIX A: PROJECTION OF THE 
EQUATIONS OF MOTION ONTO THE 
NON-ORTHOGONAL ORBITAL BASIS 



where the component Pr = Pr and by capitol P's we 
mean canonical momenta versus smah case, which wih 
mean components. To break down the LHS involves 



1. Projection onto Orbital Basis 

By projecting the equations of motion onto the orbital 
basis (n, ^, we show here that the equations of mo- 
tion depend only on the radius. 

The four equations of motion in the orbital plane are 
obtained by projecting Hamilton's equations onto the ba- 
sis vectors, as is done in celestial mechanics. For now, 
consider only the projections onto the orbital basis vec- 
tors to generate the four equations. 



rn + rn 



r • n 
p ^ 



dH 

dp 

OH 

'd^ 
dH 
dr 
OH 

dr 



n 



n 



(Al) 



To break down the LHS and RHS of the above projections 
it will be useful to write 



= (p • n)n + (fi X p) X n 
L 

= P^fi H X n 

r 

L ^ 

r 



(A2) 
(A3) 

(A4) 



p = P^n + P^A- ^r^ + . (A5) 



So, we will need n and ^, which are most directly ob- 
tained by expanding (n, ^) in the (X, Y) basis and then 
expanding (X, Y) in the Cartesian basis. So, we will need 

n = cos ^X + sin ^Y 

^ = — sin <I>X + cos $ Y 

X cos^i + sin^j 

Y = sin6>Y(— sin ^i + cos ^j) + cos ^yk 

^ = -sin^i + cos^j . (A6) 

Using 

It = ^4/ = nL^ (A7) 

^ = -^X = -I^lX (A8) 
Y = sin<9y^ = -sin i9yl]LX . (A9) 

where = Ql = SiJ/r^ from Eq. ([26|) so that we have 



+ cos $X + sin $Y 
+ Ql (^cos - sin $ sin X 



J}l_ 
sin 6y 



cos ^Y — cos ^ cos ^y k — sin <l> sin 6>y X 



(AlO) 



To take the projection of Eqs. (jAip we will also need since X lies in the equatorial plane, it is by definition 

perpendicular to J = k. From all of the above relations 
fi • X = cos ^ ^ • X = — sin ^ we obtain for use in the projections 

n-Y = sin<l> ^•Y = cos<l> x ^ ^ / a i o\ 

n • n = 

n • k = sin <l> cos ^y ^ • k = cos ^ cos ^y (All) -A /kio-/3 /kio n ^^-\A\ 

In the last step we use 

Now for <6. Taking the derivative of <l> as expressed in 
k • fi = J • (cos <I>X + sin <I>Y) = sin <I>J • Y (A12) Eq. (IMj) we have 
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-^n - sin <I>X + cos 

-<i>n + sin - cos $ sin 6>y 



sm Uy 



sin $Y + sin <l> cos 6>y k - cos <l> sin^ 6>y 



(A15) 



and using Eqs. (jAlip . we have for use in the projections 
of Eqs. dSH), 



^•n = - (^i> + J^L sin 6>y^ = - (^i> + ^]L cos 6>l) 

^•^ = . (A16) 

Now we can derive the equations of motion in the 
(r, ^, ^) coordinates. From the equations we constructed 
in section §11 C\ 



S X 



Sxp L-S. 



35i^h , (A17) 



and the projections (Eqs. ()Aip ). with ah of the the above 
vector relations we have the radial equation, 



r n = r = APr + B 
The ^ equation is found from 

r ^ = —— ^ 
op 



(A18) 
(A19) 



, L , , (S X r) • * 

r 



Look at 

(S X r) • ^ = r (S . l) . (A21) 
The <l> equation is then 

^ = A^^QlI^- cos Ol + ^ j (A22) 

where S • L is constant. Another helpful relation is 

S-L ^ ^ L L , , 

— - = J • L - - = cos6>L - - (A23) 

U (J (J 



allowing us to write the <l> equation in its final form, 

(A24) 



J 

The two conjugate momenta equations are next. We start 
with Pr'. 

p-n = i=V - -(^> + ^L cos i9l) (A25) 
r 



CPr^D^29.L 



S L 

~17 



where we have used that 

(p X S) • fi = ^ (A26) 
r 

Notice if we use Eq. (|A22p . we have 

Pr = ^CPr^D^ 31^L^ 

and last 



p-^ = P^(i> + 1^1, cos 6>L)-—r (A27) 

r J 
where we have used that 

(p X S) . ^ = S • X p) = -PrS . L (A28) 



Notice if we use Eq. (|A22p . we have a cancellation and 
(APr-r)^ = -B- = C- 

which confirms a true statement but does not provide any 
new equation of motion since we implicitly used L = 0. 
We will show in the next subection that the canonical 
momentum P^ = L and so the last equation of motion 
corresponds to i\ = 0. All four equations in the orbital 
basis are compiled in the boxed Eqs. ([28|) . 



2. Conjugate Momenta for $ and ^ 

We can show that the momentum conjugate to <l> is 
P^ = L and the momentum conjugate to ^ is = 
Lz = LcosOl' So the equations of motion (<I>,P<|>) and 
(^,P^) should be derivable from 



dH 
dP^ 
dH 
dP^ 



P^=0 
P^ = 



(A29) 



This is far more elaborate than one might guess and so 
we spend this last subsection verifing that L and Lz are 
the conjugate momenta and that the equations of motion 
derived according to ()A29p are in fact the same as those 
of Eqs. ([281). 
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We begin by showing that = L and = Lz are 
consistent with our equations before we exphcitly red- 
erive the equations of motion using (|A29p . Bear in mind 
that the variables (r, <l>, ^) and their conjugate momenta 
must be hnearly independent and so dX'^/dX^ = Sj 
where X = (r, P^, ^, ^<3E>, ^, ^ve^)- We also need to be care- 
ful to rewrite everything in terms of (r, <l>, ^) and the con- 
jugate momenta (Pr, L^Lz). Particularly, we will need to 



take the derivative of terms like cos Ol 
From Eqs. ()A29p we have that 



P^/P^ = Lz/L. 



Now 



dH _ dH dp dH Or 
dP^ ~ 'dp ' dP^ ^ ' dP^ 

. dp . dr 
P- 



dL 



dL 



(A30) 



and 



dp 

dL 

dT_ 

dL 



dn 

^dL 



Ld^ 

r dL 



(A32) 



The unit vectors fi and ^ depend on L and Lz through 
cos Ol = Lz/L. Using Eq. (|A6p and replacing sin^y with 
cosOl etc. we have 



dn 
dL 
d^ 
~dL 



. ^dY . ^cost/L ^ r\ 

sm <l>— — = - sm <l> — - — ^ - cot Ol^ 
dL L \ / 



dL 
dY 



L 

cos <l>^ = - COS ^ ^^^^^ - cot (A33) 



p — PrU + P^n 



Lr 



^ L X 



(A31) 



Using Eq. (|ATq1) for n and Eq. (IAT4|) for fi • ^, Eq. (|A3Ql) 

becomes 



^ = $ + ^ cos 6>L H — (^cot diX - ^ ^ 



Pr sin $ H cos $ j ^rn + rfi^ 



Lr - L - 

r sin <l> ( P^fi + P^fi H 4 



Taking the dot products, there are some fortuante can- 
cellations and overall we find 



<|) = <|) + ^ cos 6>L - ^ cos (9l = ^ 



(A34) 



As claimed, P^ = L is consistent. The same procedure 
for ^ with P^ = Lz yields a similarly consistent equality. 
Now, to derive the equations of motion directly from 



dH 
~dL 



dH 



PN 



dH, 



so 



dL 



dL 



(A35) 



is tedious but doable. The PN piece is straightforward if 
one first writes everything in terms of canonical variables; 
i.e., (n • p) = Pr and p^ = P^ + pj/r^ = p2 + L^/r^, 
then 



dH 



PN 



L 



dL 



(A36) 



The partial of the spin orbit piece is 
dHso _SiS d{T X p) 



dL 



6iS 



dL 
dh 

dh 



- r X 



r X 




Replace P^n with p — (L/r)^ to cancel the first cross 
product with part of the second cross product. We then 



have 



dHso 
dL 



Ol o 1 5-» 



n X 



5* 

'dL 



i dh 
^""dL 



Taking the cross products we have 



dHso ^ S_L _ ^^^^^ _ ^^^.^ 



dL 



Si 



S-L ,5i 



J -L^- cot 6*^8 



where in the last step we have used S-k = J — L^. Taking 
S • * = (J - L) • * = -L • * = -Lsin6'L gives 



dH. 



so 



dL 



S-L 5i 



cos9l {J — Lz — L cos 9l) 



= Ol — r ^ COS Ol 

-"4 ^ 



(A37) 



Added together Eqs. ()A36p and ()A37p give the equation 
of motion for <l> in Eqs. (|28|) as claimed. The 4^ equa- 
tion can be derived similarly. Since both P^ and P^ are 
constants, the Hamiltonian is cyclic in <l> and ^. 
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